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We have recently performed a relativistic 0(q 4 ) chiral expansion of the two-pion exchange AW 
potential, and here we explore its configuration space content. Interactions are determined by three 
families of diagrams, two of which involve just qa and , whereas the third one depends on empirical 
£f>l ' coefficients fixed by subthreshold nN data. In this sense, the calculation has no adjusted parameters 

and gives rise to predictions, which are tested against phenomenological potentials. The dynamical 
structure of the eight leading non-relativistic components of the interaction is investigated and, in 
most cases, found to be clearly dominated by a well defined class of diagrams. In particular, the 
central isovector and spin-orbit, spin-spin, and tensor isoscalar terms are almost completely fixed by 
just gA and /„-. The convergence of the chiral series in powers of the ratio (pion mass/nucleon mass) 
is studied as a function of the internucleon distance and, for r > 1 fm, found to be adequate for 
most components of the potential. An important exception is the dominant central isoscalar term, 
where the convergence is evident only for r > 2.5 fm. Finally, we compare the spatial behavior of 
the functions that enter the relativistic and heavy baryon formulations of the interaction and find 
that, in the region of physical interest, they differ by about 5%. 



O 1 I. INTRODUCTION 

o 

The research programme for the study of nuclear interactions was outlined more than fifty years 
f"^ , ago, in a seminal paper by Taketani, Nakamura, and Sasaki [1]. Pions, then recently detected, 

**•>. were identified as the relevant degrees of freedom for the construction of a theoretical potential. 

One pion exchanges would dominate at large distances, the exchanges of two uncorrelated pions 
would come next, and a square well could be used to simulate short range processes. It is quite 
remarkable that these ideas could stand for such a long time, survive the QCD revolution, and still 
remain as the qualitative framework of contemporary research. On the other hand, when the AW 
research programme was first established, no precise information concerning the intrinsic structure 
of pions and their interactions with nucleons was available. It took about forty years of intense 
collective work, both experimental and theoretical, for this aspect of the problem to be tamed, with 
■ the formulation of chiral perturbation theory (ChPT). 

The present day rationale for describing nuclear interactions by means of chiral symmetry is 
that low energy processes are strongly dominated by the quarks u and d and one may work with 
a two-flavor QCD. As the masses of these quarks are small in the GeV scale, one treats them as 
perturbations in a massless Lagrangian. The theory is symmetric under the Poincare group and, 
in this limit, also invariant under both isospin and chiral SU(2)x SU(2) transformations. This last 
symmetry is realized in the Nambu-Goldstone mode and the QCD vacuum contains a condensate, 
that allows collective excitations, identified as pions. The non-Abelian character of QCD prevents 
low energy perturbative calculations and, in practice, one works with chiral effective theories, in 
which point-like baryons interact by exchanging pions that have small masses. 

The one-pion exchange potential (OPEP) became definitively established in the early 1960s and 
is assumed to dominate completely AW partial waves with orbital angular momentum L > 5. Its 
mathematical form was determined in the 1950s and remains stable ever since. One has also learned 
that any ttN interaction Lagrangian, based on either pseudoscalar or pseudovector couplings, chiral 
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symmetric or not, yields the very same OPEP. Chiral symmetry is thus irrelevant for this part of 
the force, as for all single pion processes. 

The very opposite happens with the next layer of the interaction, the two-pion exchange potential 
(TPEP). This component is closely related to the 7T-/V scattering amplitude and chiral symmetry 
becomes extremely important. In the 1960s, no perturbative treatment for strong interactions 
was available [2] and potentials were constructed which incorporated nN information by means of 
dispersion relations [3]. In the same decade, chiral symmetry was being developed in a different 
framework and, with the help of current algebra techniques, low energy theorems for many pionic 
amplitudes were derived. Applications of chiral symmetry to NN interactions [4], three-body forces 
[5], and exchange currents [6] began to be performed in the 1970s. At the end of this decade, 
Weinberg [7] outlined a research programme based on the idea of ChPT. In the 1980s this theory 
was fully developed for the meson sector [8] and began to be used in the study of meson-baryon 
interactions [9]. 

The systematic use of ChPT in the study of nuclear forces began in the early 1990s, through the 
works of Weinberg [10] and Ordonez and van Kolck [11], followed by other authors [12,13]. These 
early attempts to construct a chiral TPEP considered only pion and nucleon degrees of freedom and 
gave rise to poor descriptions of NN data. Realistic potentials require other degrees of freedom, 
which were introduced in the form of deltas [14], hidden within -kN subthreshold coefficients [15,16], 
or incorporated into low energy constants (LECs) of effective Lagrangians [17-21]. In spite of 
apparent differences, there must be a rather important overlap among these various approaches. 
This is expected because the numerical values of the LECs are normally obtained from empirical 
nN subthreshold coefficients which, in turn, are largely dominated by delta intermediate states [22]. 
So, to a large extent, one is just using different languages to express the same physics. Support 
to this view comes from the fact that potentials based on deltas [14], subthreshold coefficients [23], 
or LECs [17-19,24] could produce satisfactory descriptions of asymptotic NN phase shifts, without 
free parameters. This suggests that, if one could control carefully the peculiarities of the various 
approaches, the hope of having a TPEP as unique as the OPEP could be realized. This uniqueness 
is of major theoretical importance, since it would indicate that the effective theory can indeed 
represent QCD. 

In ChPT one uses a typical scale q, set by either pion four-momenta or nucleon three-momenta, 
such that q <C 1 GeV. The leading term of the chiral TPEP is 0(q 2 ) and, at present, there are two 
independent expansions of the potential up to 0(q 4 ) in the literature. The first one is based on heavy 
baryon chiral perturbation theory (HBChPT) [17,20], where one uses non-relativistic Lagrangians 
from the very beginning and the inverse of the nucleon mass (m) as an expansion parameter. 
Relativistic corrections, needed at 0(q 4 ), are added separately [21]. The alternative calculation was 
produced recently by ourselves [25], is covariant, and results were expressed directly in terms of loop 
integrals and observable subthreshold nN coefficients. In the case of nN scattering, heavy baryon 
[26] and relativistic [27,28] results do not coincide, due to the presence of some diagrams [29,30] 
that cannot be represented by series in powers of q/m. The same class of diagrams is present in the 
TPEP and the relativistic potential also cannot be expanded in the heavy baryon series around the 
point t — 4/i 2 . If this restriction is nevertheless ignored and the q/m expansion is performed in the 
relativistic potential, one recovers most of the structure produced by the heavy baryon formalism. 
The main discrepancies take the form of both 0(q 3 ) and 0(q 4 ) terms, associated, respectively, with 
the iteration of the OPEP and the Goldberger-Treiman discrepancy. The latter could be easily 
incorporated into the heavy baryon formalism, whereas the former may derive from the particular 
definition adopted for the potential, that has to suit a given dynamical equation. In our calculation 
we treated the iterated OPEP as in Ref. [31]. 

QCD is a well defined theory and the same should happen when one works at the effective level. 
In the case of the TPEP, we consider the partial convergence between heavy baryon and covariant 
results at 0(q 4 ) as a rather welcome indication that uniqueness may not be too far ahead. The 
considerable narrowing of the theoretical discussion in the last decade represents a measure of the 
progress promoted by the systematic use of chiral symmetry, which has allowed one to understand 
the internal hierarchies of the NN potential in terms of chiral layers. Nevertheless, the question 
still remains open as to the extent this mathematical picture is backed by nature. 

The chiral picture may be assessed by comparing theoretical and empirical phase shifts. In the 
case of peripheral waves, this can be done perturbatively, without resorting to free parameters, and 
the main trends are well reproduced [23,17,24]. However, these waves are small, error bars are 
important, and the test is not very stringent. In order to include inner waves, which are larger, 
one has to use dynamical equations, but this demands the regularization of the chiral potential by 
means of form factors, since it diverges at large momenta. The problem with form factors is that 
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their use amounts to extending the dynamical content of the potential to higher orders in q/m by 
means of free parameters, in a way which is not controlled by chiral symmetry. This blurs the chiral 
content of the potential and success in reproducing data, welcome as it is, also does not represent a 
very stringent test. 

As an interesting alternative, we may assume that the chiral potential, calculated at a given 
order, determines completely the interaction from a radius R onwards and then use it as an input 
in phase shift analyses. This would just amount to extending to the TPEP a procedure which has 
already been used for a long time in the OPEP. For the latter, this idea has proved to be reliable in 
the elastic regime and for waves with L > 5. From the standpoint of the symmetry, this happens 
because chiral corrections are short ranged and one sees just the leading contribution through this 
window, irrespectively of the order in q one is working at. In the case of the TPEP, the corresponding 
problem is much more complex and not fully understood. Works along this line have already been 
performed by the Nijmegen group [18], who claim that a 0(q 3 ) potential is effective for distances 
smaller than 2 fm. However, its conclusions are disputed by Entem and Machleidt [32] and the 
situation remains unclear. 

The present paper is motivated by the feeling that the quantitative aspects of chiral hierarchies 
need to be clearly understood if the TPEP is ever to become a reliable tool to be used in phase shift 
analyses. Our study is based on the configuration space version of the 0(q 4 ) potential produced 
in Ref. [25] and organized as follows. In Sec. II we discuss the dynamical content of the TPEP, 
which is given by a set of Feynman diagrams, organized into three families. The explicit expressions 
and corresponding figures for the various components of the potential are given in Sec. III. As the 
way chiral symmetry is implemented varies with the family considered, in Sec. IV we discuss how 
dynamics is mapped into the final form of the potential and show that the importance of the LECs 
is rather channel dependent. Sec. V deals with the convergence of the chiral series and in Sec. VI we 
discuss the main differences between the relativistic and heavy baryon approaches to the potential. 
Finally, conclusions are presented in Sec. VII. 



II. DYNAMICS 




FIG. 1. Two pion exchange amplitude. 

The dynamical content of the relativistic 0(q 4 ) chiral TPEP is determined by the effective La- 
grangian 

£ e// =4 2) +4J ) +^ ) +4 v 3 \ (2-1) 
where and cffl describe pion-pion and pion-nucleon interactions at 0(q n ). Other degrees of 

(2) 

freedom are implicitly taken into account by means of the LECs d and di, present in C N and 
■ The use of covariant Feynman rules with vertices derived from this Lagrangian allows the 
construction of the T matrix T$ , which describes the on-shell process N(pi) N(p2) — > N(p' 1 ) N(p' 2 ) 
and contains two intermediate pions, as represented in Fig. 1. The potential is obtained by going 
to the center of mass frame and subtracting the iterated OPEP, in order to avoid double counting 
when it is used in the Lippmann-Schwinger equation. 
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FIG. 2. Dynamical structure of the TPEP. The first two diagrams of family I correspond to the products of Born ttN 
amplitudes, the third one represents the iteration of the OPEP and the next three involve contact interactions associated 
with the Weinberg- Tomozawa vertex. The diagrams of family II describe medium range effects due to pion-pion correlations. 
Interactions represented by family III are triangles and bubbles, involving nN subthreshold coefficients, indicated by the large 
black dots. 

The dynamical content of the amplitude T^p is given by the diagrams of Fig. 2. Their full 
evaluation produces amplitudes containing many different loop integrals, which are interconnected. 
The chiral orders of the potential are extracted by exploring as much as possible the mathematical 
relations among the various loop integrals. As the use of these results represents an important step 
in the determination of the potential, in the Appendix we display their accuracy in configuration 
space. 

The processes given in Fig. 2 are organized into three different families. The first one corresponds 
to the minimal realization of chiral symmetry [13], includes the subtraction of the iterated OPEP 
and involves only the pion-nucleon interactions given by , with the constants m, g and / w used at 



their physical values. The second family contains two-pion correlations in the t channel, determined 

,(i) 

-N 



by and . Finally, the last family includes chiral corrections representing either higher order 



processes or other degrees of freedom, hidden into the LECs of jC^ and jC^ . 

This theoretical structure has been fully incorporated into our recent evaluation of the amplitude 
T^p , Ref. [25]. In that work we have performed a two-step calculation, using the fact that the 
NN interaction is closely associated with the off-shell nN amplitude. This allows one to use many 
of the results derived by Becher and Leutwyler [28] (BL) for the nN amplitude as inputs into the 
evaluation of the NN potential. Moreover, it clarifies the relationship between the chiral orders of 
the AW and nN amplitudes. Using Fig. 1, we write the 0{q n ) expansion of Ttp as 

r (n) _ _ j_ 1 / I d_Q I \p r T (i) ,(1) r T (m),(2) ,„ „x 

TP ~ 2! (27r) 4 J 1 [k 2 — fj, 2 ] [k 12 — y?\ J 2^ l J *"J , W 

where [T^]^ is the nN amplitude for nucleon (i) expanded at order 0(q m ). The factor within 
curly brackets in the integrand is O(q ) whereas the leading term in T^n, as given by the Weinberg- 
Tomozawa theorem [33,34], is O(q). Thus T^p requires T n N up to 0(q n ^ 1 ). 

This result is important regarding the numerical values of the LECs to be used in the determina- 
tion of the TPEP, which depend on the chiral order one is working at [35]. These constants are not 
observables and must be obtained from empirical quantities such as, for instance, nN subthreshold 
coefficients. In the case of our 0(q 4 ) TPEP, consistency demands the use of LECs determined from 
T^N at 0(q 3 ). 

Finally, a further motivation for deriving the TPEP from the intermediate nN amplitude is that 
this stresses the continuity of present developments with the seminal works of the Paris [3] and 
Stony Brook [4] groups, produced more than three decades ago. For this very reason, one becomes 
better prepared to understand the specific role played by ChPT in this problem. 
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III. CONFIGURATION SPACE POTENTIAL 



The configuration space Schrodinger equation is a rather useful tool for calculating low energy 
nuclear processes. In principle, the 0(q 4 ) r-space potential could be obtained by just performing 
the Fourier transform of our center of mass p-space potential, which is written as 1 



t C m = 3t + +2r (1) -r (2) t 



(3.1) 



with 



LS ,± 



T ,± . ^SS ,± 



(3.2) 



and n LS = i (o- (1) +o- (2) ) ■ q X z/4, fi T = -q 2 (3o-« • qa™ • «f-<r« ■ o-< 2 >), Sl ss = q 2 er« • a {2 \ 
and J1q = cr' 1 ^ ■ qxz a 1 ' 2 ' 1 ■ qxz. However, this leads to expressions that contain non-local terms, 
due to the energy dependence of the profile functions U. In order to avoid this kind of complication, 
we expand the potential in the non-local operators and keep only local and spin-orbit contributions. 
In this approximation, the configuration space potential becomes 

V(r) = (v£ + Vi s Mls + V£ Q,t + V^s ^ss) 

+ r (1) • r (2) (ye + V£ g n LS + Vt n T + Vgs n ss ) 

with f! LS =I-( ff W+(7( 2 »)/2, fi T = 3o-« r ■ r - -a™, Q S s = <t W ■ t {2) ■ 



The radial functions are given by 



V£(r) = T ± U${x), 



Vts(r)=r i 



fj 2 1 d 

m 2 x dx 



(3.3) 

(3.4) 
(3.5) 



dx 2 x dx 



U±(x), 



Vssir) = -r* ^ 



dx 2 x dx 



(3.6) 
(3.7) 



where r + = 3, r = 2, x = fir, and 



Vf{x) = - J 



k ik x *± 



(2tt)3 



tf(fc), 



J = {C,LS, T, SS} 



with fe = q/fi. This allows the potential to be expressed in terms of dimensionless configuration 
space Feynman integrals, denoted by S, and related to the functions II of Ref. [25] by 



S(x) 



I 



3fc e ikx U(k). 



(2k) 



(3.8) 



(3.9) 



Using the results of Sec. IX of Ref. [25], we have the expansions 2 



tj+ _ M 3 m 2 
c 256tt 2 /* 



{g\ (1 + 4 A GT ) (1- V 2 /2) 2 (S x - S b ) 



flA 



(1- V 2 /2) [-g\ (2S a + V 2 S t ) + 8 (6+ + 8+ V 2 ) & 



+ 



■^f! (1-2V 2 ) (-4,7(1- V 2 /2)5 t + (1- V 2 /2) 2 S tt ) + ^V 2 (S X + S b ) 

1D7T J-7r 4 



J In this result, the (+) and (-) upper labels indicate, respectively, terms arising from the isospin even and odd irN 
subamplitudes. 

2 In writing these expressions, we did not consider the relativistic normalization factor, proportional to m/E. 



5 



[slV 4 - 4g 2 A ((5+ + 8+ V 2 )V 2 + 4.(1 - 2V 2 /3 + V 4 /6)) 
+8 (5++S+V 2 + (5+/3)(l - V 2 /4)) 2 + §(<5! + o) 2 (l-V 2 ) 2 ] S e } , 



U7.a = 



q 2 2 

fi m g A 



9 2 a 



{g 2 A [(1 - V 2 /2)(5 b - S t ) - (3/2 - 5V 2 /8) S a ] 

16 



: ' " ~ 128tt 2 / 4 

(1 + 2V 2 - V 4 /2) (S* x + ft) + (2 <?iV 



" j4,(l"V a /4))ft 



(3.10) 



(3.11) 



U. 



"^f = -S5f {-Pi(l + 4AGT)(1-V 2 /4)[S X +S 6 ] 



^ [(1 - V 2 /2)(ft - S 6 ) + (1 - V 2 /4)ft] 



^(l-V 2 /2) 2 S x +^+(l-V 2 /4)ft 



and 



TJ- - M 3 m 2 

° 16^ 2 / 4 



|| (1 + 4 A GT ) (1 - V 2 /2) 2 (S X + ft) 



-J [si(l + 4A GT )-si(l + 2A GT )] (1- V 2 /2)ft 

+ ^ [ffl (1 + 4 A GT ) - 2g\ (1 + 2 A GT ) + l] (1 - V 2 /4)ft 

^(l - V 2 /2) [si(ft - V 2 ft/2) + ( S i - 1)(1 - V 2 /2)ft] 



2 f.2 



^(1 - V 2 /2) [-( S i - l)V 2 /8 + 5 - + 8 m V 2 + 5 W (1 - V 2 /4)/3 + A»V 2 /4] 



(g ^ 6 1} (1 - V 2 /4) [J - + 5 M V 2 + 3^,(1 - V 2 /4)/5 + ^ooV 2 /4] J ft 

[(2si(l - 5V 2 /6 + V 4 /5) + 4(si - 1) 2 (1 - 3V 2 /8 + V 4 /32)/9 



64tt2/J 



-4si(fli - 1)(1 - 29V 2 /72 + 7V 4 /144)) 5, 

(sl(l - V 2 /2) 2 - 2si(si - 1)(1 - V 2 /2)(l - V 2 /4)/3 + (si - 1) 2 (1 - V 2 /4) 2 ) S u ] 



2 f|(l-V 2 /2) 2 (V 2 /4)(5 x -ft) 



U7z = - 



fi 3 m 2 
128tt 2 /| 



{si [(3/2 - 5V 2 /8)ft - (1 - V 2 /2)(ft + ft)] 



+2si(si - l)(l - V 2 /4)ft 



[(si - 1) 2 (1 - V 2 /4)/2 + 4si Poo (1 - V 2 /2) - 4 (si - 1)^,(1 - V 2 /4)/3] ft 



^■(1- V 2 /2) 2 (ft -ft) 



(3.12) 



(3.13) 
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2 4 

m g A 



8^71 h 27r d-v 



+ (1 - V 2 /4) 2 S it 



}• 



(3.14) 



-1g\{g\-\-2 



1536tt 2 /* 

i OT )(l-V 2 /4)S 



{<?![(! -V 2 /2)S b + (l-V 2 /4)S a ] 



+ 



- 1 - 2/3 M )(l - V 2 /2) - (g\ - 1 - 2 Ak.) 2 (1 - V 2 /4)/3l & 



|^[-2^ ( l-V 2 /4)S t + (l-V 2 /4) 2 S tt ] 



where the Laplacians act on the variable x. The chiral orders of the various radial functions may be 
read directly from the combination of Eqs. (3.4)-(3.7) and (3.10)-(3.15). Their relative importances 
will be discussed in detail in Sec. V. We have expressed our results in terms of the axial coupling 
gA- If one wants, they may be rewritten using the ttN coupling constant g, by means of the relation 
g = (1 + Agt) gAm/fjr, where Act is the so called Goldberger-Treiman discrepancy. 

The parameters Sfj and entering these expressions are determined by subthreshold ttN coeffi- 
cients or, alternatively, by the LECs of the effective Lagrangian, according to the results presented 
in Sec. V of Ref. [25]. Their empirical values are reproduced in Table I. 



(3.15) 



TABLE I. Dimensionless subthreshold coefficients; definitions are the same as in Ref. [25]. 



*0 + 








-4.72 


3.34 


4.15 


-10.57 








fa 


7.02 


-3.35 


-2.05 


5.04 
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The eight functions £» which carry the spatial dependence of the potential are dimensionless and 
given by 



Si 

S a 



Ki(2x) 



-2x 



2x 2 ' 

1 Pi 

da / db 

o Jo 



(1 - b) 2m/ 'n e 



-S t x 



A? 



4-kx ' 



A 2 = a(l-a)(l-bf , 



61 = [(1 -b) + b 2 m 2 /^} /A? , 

fdb fta^Ljglf ^ 
o Jo x 



8tt0x ' 



A" =o(l-a)(l-6) , 

6> 2 = [(1 - ab) + a 2 b 2 m 2 /^ 2 ] /A 2 , 

~ 2 b 4to 2 /a« 2 e~ e " x 



poo pi pi ^2 

Sb = I dc db da — 
Ji Jo Jo 



8iv8 b ' 



A 2 = o(l - a)(l - 6) - aV(l - c 2 )/4 , 
6» 2 = [(1 - ab) + a 2 b 2 c 2 (m/n) 2 ] /A 2 , 
1 f 1 « 3 & 2 4m 3 / M 3 



/■OO /■! PL ^ 

Sb = I dc I db I da — 

Ji Jo Jo 



Su = 



1 db^ {1 ~ a){1 ~ b) 1 1 



da / d& 



(6 + a) z(I-I) 
G(V/a)G(V/6) 1 



^ v 2 e -2x/a / 4 \ 2 e ~ 2x / b 



WJ _ r VP" J 



(4m^) 2 

' ; ""; a 2 b 2 (« + &) 

where K\ (x) is the modified Bessel function and 



1 1 r 



-2x/ a 



-2x/b 



G(t') 



V) 



arctan 



y/(4m 2 -t')(t' - 4^ 2 ) 



In Figs. 3(a)-3(h) we display the numerical predictions of our TPEP (full line), obtained by using 
the parameters 5^ and (if- given in table I, fixed by the nN subthreshold coefficients of Ref. [22]. 
As we will discuss in the sequence, our chiral TPEP is theoretically reliable for large distances and 
definitely not valid for internucleon separations smaller than 1 fm (shaded area). For the sake of 
producing a feeling for the phenomenological implications of these results, we also plot the medium 
range components of the Avl4 [36] and Avl8 [37] versions of the Argonne potential (dotted and 
dashed lines respectively). 

The central isoscalar component of the nuclear force is by far the most important one and the 
fact that the chiral prediction is consistent with both Argonne versions is rather reassuring. 3 The 
assessment of the other components is more difficult, since there are important variations between 
the Avl4 and Avl8 results. In the cases of Vg S , Vf , V^, where these variations do not involve 
signs, it is possible to note a qualitative agreement with the behavior of the chiral TPEP. The curves 
for Vg~ s , Vq and V^T are not far from those of Avl8 whereas coincides with the Avl4 prediction. 

In order to complete the long-distance description of the NN potential, one has to include the 
OPEP, which contributes only to V^T and Vg S , through the following expressions: 



(3.16) 
(3.17) 
(3.18) 



(3.19) 



(3.20) 



(3.21) 
(3.22) 
(3.23) 

(3.24) 



3 The leading structure of V<4~ was discussed in Ref. [38]. 
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V t \ OPEP + 2 &gt) (, 2 + 3, + 3) V - 0-25) 



\opep e 48nm 2 



x- 1 



These components, which dominate at large distances, are shown in Figs. 4. (a) and 4.(b), together 
with the corresponding TPEP contributions. The influence of the TPEP only becomes significant 
in Vf for r < 2 fm, and in Vg S , for r < 3 fm. 
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FIG. 3(a). Central isoescalar component. 
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FIG. 3(b). Spin-orbit isoescalar component. 
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FIG. 3(c). Tensor isoescalar component. 




FIG. 3(d). Spin-spin isoescalar component. 
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FIG. 3(e). Central isovector component. 
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FIG. 3(f). Spin-orbit isovector component. 
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FIG. 3(g). Tensor isovector component. 
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FIG. 3(h). Spin-spin isovector component. 
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FIG. 4. (a). OPEP and TPEP contributions to the tensor isovector potential. 




IV. INTERNAL DYNAMICS 



In this section we discuss the relative importance of the contributions originating from the three 
families of diagrams presented in Fig. 2. This is motivated by the fact that the chiral description of 
the TPEP consists of a well defined field theoretical structure which depends on external parameters 
representing masses (/i, m), coupling constants (f n , gA), and LECs (cj, In order to be able to 
obtain predictions, one has to feed the mathematical structure with the empirical values of these 
parameters. 

The constants present in the 0(q 4 ) potential may be divided into two classes, according to their 
numerical accuracy. The values of fx, m, f v , and gA entering and may be considered as 
being very precise for the purposes of determining the TPEP. On the other hand, the constants a 
and di that appear in and jCjy need to be extracted from ttN subthreshold coefficients by means 
of dispersion relations and hence may contain both experimental and theoretical uncertainties. This 
means that, in the case of the interactions given in Fig. 2, predictions from families I and II are 
very reliable whereas those associated with family III may be less so. For this reason it is important 
to establish how the results discussed in the preceding section depend on the various families of 
diagrams. 

In order to assess the importance of each family we show, in Fig. 5, their relative contributions 
to the components of the TPEP. A general pattern one can observe is that two-loop contributions 
(family II) are negligible and, in particular, exactly zero for V^ s , Vg S , and ■ The various profile 
functions are neatly dominated by either family I or III. The former, which is very precise, dominates 
the channels Vq , V^ s , Vg~ s , and and a modification on the values of the LECs would hardly 
influence the corresponding curves. This can be viewed as a strong constraint on the construction of 
phenomenological potentials. For the remaining channels, this condition is somehow relaxed, since 
they are dominated by the diagrams of family III. If one wishes, the freedom in these channels may 
be used to fix experimentally the LECs by means of NN data. 
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Fig.5(a) V c + (partial) / V c + (total) 



Fig.5(b) V L + S (partial) / V L + S (total) 
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Fig.5(c) V T + (partial) / V T + (total) 



Fig.5(d) V s + S (partial) /V s + S (total) 
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Fig.5(e) V c ~ (partial) / V c ~ (total) Fig.5(f) V L " S (partial) / V L " S (total) 



2 








10 


/ 

/ 

- / 






1 










- / 

- / 
1 

/ 

/ 












5 





















1 










-5 


- 


/ 

/ 

1 

1 

1 

1 






family I 

family II 

family III 






family 1 

family II 

family III 




2 








-10 




i / 
i . / 





2468 10 2468 10 

r (fm) r (fm) 



Fig.5(g) V T ~ (partial) / V T ~ (total) Fig.5(h) V s " s (partial) / V s " s (total) 




2468 10 2468 10 

r (fm) r (fm) 

FIG. 5. Relative weight of each family in the TPEP, obtained by dividing the partial contributions by the full result. 

V. CHIRAL STRUCTURE 

In this section, we discuss chiral scales. In the case of the central components, these scales can be 
read directly from the functions Uq, given by Eqs. (3.10) and (3.13). For the other terms, there is a 
factor (/i 2 /m 2 ) in the relation between V^ 1 and U ± , arising from the non-relativistic expansion of the 
Dirac spinors, (3.5)-(3.7). Thus, in the 0(q 4 ) potential, one expands the corresponding functions 
U ± up to 0(<? 2 ). 

The leading term of the chiral TPEP is 0(q 2 ) and our results are written as sums of 0(q 2 ), 0(q :i ), 
and 0(q 4 ) contributions. In the cases of VjJ", Vg~ s , and Vq , this structure is mapped directly into 
the corresponding profile functions. The other components begin at 0(q 3 ). 

In p-space, the chiral series involves nucleon three-momenta, assumed to be small. This means 
that, in r-space, the chiral structure should become apparent at large distances. In order to check 
this, in Fig. 6 we show the ratios of the chiral layers for the various components of the potential. In 
all figures it is possible to note, at large distances, a rather well defined chiral hierarchy. Corrections 
are always smaller than the terms they correct. On the other hand, this hierarchy tends to break 
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down when distances decrease and we assume that our results are not physical for r < 1 fm. In two 
cases, namely, Vq and V£ s , corrections are large within the region of physical interest. 
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FIG. 6. Relative contribution of each chiral order to the TPEP. The point in the curve where the ratio is 0.5 is indicated by 
a black dot, for the sake of guiding the eye. 



VI. THE HEAVY BARYON APPROXIMATION 



The relativistic potential is expressed by Eqs. (3.4)-(3.7), (3.10)-(3.15), and involves eight basic 
functions, denoted generically by S(x). They are given by Eqs. (3.16)-(3.23), and represent bubble, 
triangle, crossed box, planar box, double bubble, and double triangle diagrams. These functions 
have been derived by means of covariant techniques and correspond to the signature of relativity 
in this problem. Only the bubble integral can be evaluated analytically and the other ones are 
not homogeneous functions of either the pion mass or external three-momenta. In general, the 
expansion of the functions S(x) in powers of q/m is not mathematically defined. However, as 
discussed by Ellis and Tang [29], if one forces such an expansion, one recovers formally the results 
of HBChPT. In Ref. [25] we have expanded our 0(q 4 ) p-space relativistic potential in this way and 
obtained (inequivalent) expressions that reproduce most of the standard 0(q 4 ) HBChPT results 
[17,20,21,24]. Differences are due to the Goldberger-Treiman discrepancy and to the procedure 
adopted for subtracting the iterated OPEP. In this section, we discuss the numerical implications 
of the heavy baryon approximation in configuration space. 

We begin by considering the triangle integral St, given by Eq. (3.18), that can also be expressed 
as [27] 



Im 7 (t') 



with 



Im y(t') 



arctan 



yj (Am? - t')(t' - Afi 2 



87r^'(4m 2 -t') 
The heavy baryon approximation consists in writing 



V - 2^ 2 



Im 7 (t') 



1 



167rm\/F 



arctan 



2m y/t' - 4fi 2 
V - 2^ 2 



and treating formally the argument 



2m J V - 4fi 2 
ex. — 

V - 2 M 2 

as being 0(g -1 ). This would suggest that one could use the result arctan a = n/2 — l/a + l/3a 3 + ■ ■ ■ 
in order to derive the heavy baryon expansion of the triangle integral. Recently BL [27] have 
discussed the properties of the spectral representation based on Eq. (6.4) and remarked that the 
series for arctan which underlies the heavy baryon approximation is valid only in the domain |a| > 1. 

For | a | < 1 one should use arctan a = a — a 3 /3H , but this corresponds to an expansion in inverse 

powers of q. They showed 4 that a suitable representation for St is 



St 



8m 
H 



°— dt' 



4,^ 



i 



167rm\/F 



arctan 



2m y / t' - 4^2 
V - 2^2 



(6.1) 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



n _ (t' - 2^ 2 ) 
2 2m v /t'-4^ 2 



arctan ■ 



ly/t'-A/X 2 



(6.6) 



They worked in momentum space. 
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The heavy baryon approximation consists in keeping only the first bracket in the integrand. How- 
ever, this does not cover the region t' ~ 4/i 2 , where the second term dominates. As a consequence, 
the heavy baryon approximation of St, which reads 



sr = 



2x 2 



+ 



^--^5 [xKo(2x) + Ki(2x)]+ 
2m nx z 



NLO 



(6.7) 



is not suitable for all values of x, as observed numerically in our previous work [38]. The exponential 
in the integrand of Eq. (6.1) shows clearly that, for large values of x, results are dominated by the 
lower end of the integration. Thus, a good description of St at large distances requires a decent 
representation for Im "lit') near t' = 4fi 2 . 
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FIG. 7. The heavy baryon expansion of the triangle integral, given by Eq. (6.7), and the relativistic BL correction (St ), 
divided by Sf L . 



In Fig. 7 we display the ratios of the various terms of Eq. (6.6) by Eq. (6.5). Inspecting this 
figure one learns that the first two terms of the heavy baryon series do not represent well the full 
result. In order to have a good description of Sf L at large distances one has to add S\ h which, as 
pointed out by BL, cannot be obtained through the heavy baryon series. 

An advantage of the heavy baryon formalism is that it gives rise to power counting, which is 
absent in relativistic baryon ChPT based on dimensional regularization [9]. In order to overcome 
this difficulty, BL proposed a new regularization scheme, based on a previous work by Ellis and 
Tang [29] . The so called Infrared Regularization (IR) respects the correct analytic structure around 
the point t' = 4/i 2 , is manifestly Lorentz invariant, and gives rise to power counting. 

In the case of the triangle integral, the infrared regularized expressions reads 



da 



db 



(1 - 6) 2m/ ix 
A? 



Atvx 



with A 2 and 6 t given by Eq. (3.18). 

In Fig. 8 we compare the infrared regularized triangle integral (St) with that given by Eq. 
(3.18), obtained through dimensional regularization. For comparison, we also plot the results of 
the heavy baryon formulation at LO and NLO. The relativistic versions of the triangle integral 
are numerically identical for r > 1.5 fm, indicating that the form of the regularization procedure is 
irrelevant in the region of physical interest. 



(6.8) 
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FIG. 8. The heavy baryon approximation of the triangle integral. 

The discussion about the triangle integral may be extended to the functions S x , St, and Sb, 
associated with crossed box and planar box diagrams, and given by Eqs. (3.19)-(3.21). The heavy 
baryon approximations for these results read [25] 



Si 



Q HB 



aHB ^ e 



4.r 2 



K (2x)- 



K (2x) 





e~ 2 * 




.m_ 


2x 
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J± 
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.m_ 


4x 


.m. 



{^ij [2xK {2x) + (1-x 2 )A' 1 (2 : e)]}h 
[2xK (2x) + {1-x 2 )Ki{2x)]} 



1J- 

I 3tvx 2 



+ ■ 



2 'xK (2x) + K 1 (2x)]+ [ V ' 



Sttx 2 



32x 



(6.9) 
(6.10) 
(6.11) 



The quality of these heavy baryon approximations may be assessed in Figs. 9-11, where the 
partial sums are divided by the relativistic result, Eqs. (3.19)-(3.21). 
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FIG. 9. The heavy baryon approximation of S x - The partial sums are divided by the relativistic result, Eq. (3.19). 
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FIG. 10. The heavy baryon approximation of Sb- The partial sums are divided by the relativistic result, Eq. (3.20). 
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FIG. 11. The heavy baryon approximation of Sb- The partial sums are divided by the relativistic result, Eq. (3.21). 



VII. CONCLUSIONS 

In this work we have studied the main features, in configuration space, of a relativistic 0(q 4 ) 
expansion of the two-pion exchange nucleon-nucleon potential derived recently by ourselves [25]. 
Chiral symmetry provides a mathematical structure for the potential, that has to be fed with 
numerical values for fi, m, / T , gA, and the LECs d and di. The main source of uncertainty are the 
values of those LECs, which need to be extracted from ttN scattering data. 

The profile functions for the various non relativistic components of the potential were compared 
with two phenomenological versions produced by the Argonne group. One finds good agreement 
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with their central scalar term, which dominates the NN interaction. In all cases in which the signs 
of the Argonne potentials coincide, there is a qualitative agreement with our results. 

In order to check how empirical uncertainties can affect numerical results, we have studied the 
dynamical content of the potential in terms of families of diagrams associated with either the [C^ + 
C^] or [Cffl + C$] pieces of the effective Lagrangian. In all but one cases, dynamics is clearly 
dominated by one of these interactions. In particular, V^ s , , V<ts, and Vq are dominated by 
[C^ + C$] and hence fixed by the values of qa and only. The components VJ4", , and Vg S , 
on the other hand, are dominated by [£}/ + Cf/} and their numerical values may be affected by the 
less certain LECs c; and dj. 

Most components of the potential are given as sums of 0(q 2 ), 0{q 3 ), and 0(q 4 ) terms. The relative 
weights of these terms of the chiral series have been investigated and one finds good convergence 
at large distances. However, there are two cases, namely, V@ and V^ s , where convergence is not 
evident in the region of physical interest. We intend to deal with this problem elsewhere. 

Finally, the relationship between relativistic and heavy baryon results has been discussed. On 
the purely conceptual side, the view seems to be well accepted nowadays that they cannot be fully 
equivalent. This is indeed the case and the numerical implications of this statement in configuration 
space were found to be of the order of 5%. 
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APPENDIX 

Several chiral calculations of the TPEP were produced in the last decade. As we pointed out 
in the introduction we expect, in the spirit of effective theories of QCD, that all these calculations 
should eventually converge to a single result. 

It is in this conceptual framework that we discuss here the relationship between the present work 
and its earlier versions, published between 1994 and 1997 [13,16]. In Ref. [25], the 0(q 4 ) expansion 
of the TPEP was performed in three steps. In step 1, we derived full amplitudes, by using standard 
covariant techniques, to evaluate the diagrams of Fig. 2. At this stage, results were quite similar 
to those of Ref. [16], although not identical, as we discuss in the sequence. A handicap of the 
full amplitudes is that they involve several cancellations and do not exhibit chiral scales explicitly. 
Therefore, in step 2 we derived intermediate results, that show these scales, by just rewriting the full 
amplitudes with the help of exact relations among Feynman integrals. We subsequently neglected 
short distance terms and, in that region, full and intermediate results became no longer identical. 
The transmogrification of the potential was based on the following relations 5 : 

Relation 01:S<° 00) = I (l - V 2 /4) Sl° c 00) +••■ , 
Relation 02:3^ = ^ (l - V 2 /4) 2 S<° 00) + • • • , 

Relation 03: [l - (^/m) 2 V 2 /4] S<° 01) = S£°°> - £ (l - V 2 / 2 ) S^ 00 ' + • • • , 
Relation 04: [l - ( M /m) 2 V 2 /4] S<°° 2) + Si° c 00) = - ^ (l - V 2 /2) Si° c 01) +■■■, 
Relation 05:5^ = \ (l - V 2 /4) Si° c 00) + £ (l - V 2 / 2 ) 
Relation 06: [l - ( M /m) 2 V 2 /4] ST } = SlT - (l - V 2 / 2 ) ST } +■■■, 



5 For the complete details of the notation, please see Ref. [25]. 
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Relation 07:S£ 00) = -S<° 01) + ■ ■ ■ , 

Relation 08: [l - ( M /m) 2 V 2 /4] S<° 02) + 5<f 0) = V 2 /2) ST } 

Relation 09: [l - ( M /m) 2 V 2 /4] 2 + [l - W™) 2 V 2 /4] Si° s 00) = A (l " V 2 /2) 2 ^ 00) 



+ [1 - ( M /m) 2 V 2 /4] S<° 01) - ^ (1 - V 2 /2) S^ ^ + ■ 

Relation 10:S<° 00) = (l - V 2 /4) S<° 00) + (l - V 2 /2) Si° s 01) + ■ 

1 



Relation U-.-Sg** = + (l - V 2 /2) S, 
2m v 



(000) 
reg 



y/1- Wm)2V 2 /4 
Relation 12: [l - ( M /m) 2 V 2 /4] S<° e ° 9 2) + S<° e ° 9 0) = S<° 01) + • • • , 
Relation 13:^° e ° 9 0) = S^ r > + JL (l - V 2 /2) S<™> + • • • , 



S a + -- 



Relation 14: 



2m 

-/^(l-V 2 /2) 2 5 r 9 o) 

4m/ v 



£(l-VV2)^-^(l-V 2 /2) 



Relation 15: S, 



(000) 



-(1 



V 2 /4)^r-^(l-V72)S res 



y/1 - (p/mY V 2 /4 

(010) , 



^(000) + . . 



In these expressions, the ellipses indicate short range terms, which have been neglected. In order 
to produce a feeling for the accuracy of these approximations, in table II we display the quantity 
Aj(r) = 1 — Ri(r) / Li(r)\, where Li(r) Ri(r) are respectively the values of the left and right hand 
sides of relation i at point r. Inspection of this table shows that, although discrepancies may be 
large at short distances, in all cases they remain below 1% for r > 1.5 fm. 

TABLE II. Approximations made in relations among integrals. 



(A7) 
(A8) 

(A9) 
(A10) 
(All) 

(A12) 
(A13) 

(A14) 
(A15) 





r = 0.1 fm 


r — 0.5 fm 


r = 1.0 fm 


r = 1.5 fm 


Rel. 01 


0.000019 


0.000000 


0.000000 


0.000000 


Rel. 02 


0.000004 


0.000000 


0.000000 


0.000000 


Rel. 03 


0.004272 


0.000001 


0.000000 


0.000000 


Rel. 04 


0.000433 


0.000000 


0.000000 


0.000000 


Rel. 05 


0.002478 


0.000001 


0.000000 


0.000000 


Rel. 06 


0.056091 


0.000002 


0.000000 


0.000000 


Rel. 07 


0.000000 


0.000000 


0.000000 


0.000000 


Rel. 08 


0.005502 


0.000010 


0.000000 


0.000000 


Rel. 09 


0.271267 


0.000347 


0.000000 


0.000000 


Rel. 10 


0.034881 


0.000009 


0.000000 


0.000000 


Rel. 11 


0.958301 


0.096190 


0.000076 


0.000312 


Rel. 12 


0.000096 


0.000030 


0.000001 


0.000000 


Rel. 13 


0.855452 


0.134455 


0.006010 


0.000168 


Rel. 14 


1.007437 


1.105614 


0.066393 


0.006432 


Rel. 15 


1.554767 


0.417655 


0.018483 


0.000513 
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Finally, in step 3, we obtain the 0(q 4 ) expansion of the TPEP, given in section III, by truncating 
the results of step 2 at that order. 

We now compare the results of this work with those from earlier versions. Our 1994 paper [13] 
dealt with the evaluation of the diagrams given in family I of Fig. 2, which corresponds to the 
minimal realization of chiral symmetry in the TPEP. The main differences with our present results 
concern second order corrections, due to the way variable W , representing the total CM energy, was 
approximated in the planar box diagram. In 1994 paper we used W = 2m, following Partovi and 
Lomon [31]. We no longer perform this crude approximation. 

In our 1997 paper [16] we have calculated the diagrams shown in family III of Fig. 2 and results can 
be directly related with those of the present work, provided one establishes the connection between 
the two notations. For instance, the central isoscalar potential was formerly written as 

^Ij/J = _ 4^2 {^m 01 ™™ [^ S ( 2m >") ^^T(2m,n)] + ^ti a rnnS B(2k+2m,i+n) 
+ g 2 ^-Pmn [25 , fl( 2m + l in ) + 2S^(2m+l 1 n)] + a ktPmnS^B(2k + 2m+l ,l+n) 



+ PktfimnSB{2k+2m+2,t+n) + 9 m Pmn ^^T(2m+l,n) + PkePmnS%^ 2 k+2m+2,e+n) J ' (A16) 

where S are Feynman integrals from Ref. [16] and atm and j3tm are linear combinations of 7vN 
subthreshold coefficients. In order to recast the old results in the form adopted in this work, one 
may use the relations 

5 fl( o,„) = ^(l-V 2 /4) n ST (A17) 

S B( 2, n) = S| (0in) = Sl (hn) = t^L (i _ v 2 /4)" ST (A18) 

Sb ( 4,„) = 3S| (2 , n) = S£ (3 ,„) = (3/2)S^ n) = t^L (i _ v 2 /4) " C° (A19) 

S^,o, n) = ^(l-V 2 /4)" (-)ST (A20) 

S T(2,n) = SY(l,n) + S T(l,n) (A21) 

S%i,n) = ^ (1 - V 2 /4)" (-) + 2ST] (A22) 

^ (1 ,„)=-^(l-V 2 /4)" (-)ST (A23) 
The parameters a^ n and /3^ n are related to subthreshold coefficients by 

<4> = A 4 (4o + V<Si + 16/40 ( A24 ) 

ati = 4/ (4i + 8A* 2 4 2 ) (A25) 

a+ 2 = 16a* 5 (4 2 ) (A26) 

a+o = A* 3 [4o - &S, + 4/ (4 - &3i )] ( A27) 

a+ = 4 M B (d+ - 6+) (A28) 

/?oo = A* 3 (bto + 4/Ali + 16aA! 2 ) (A29) 

/?oi= 4^(^1+8^) (A30) 

13+ = 16a* 7 (& ( | 2 ) (A31) 
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Just as an example, using these rules in the case of the triangle contribution to Eq. (A16) and 
truncating at 0(q 4 ), we find 



Vc \iu - 4 ^ 2 5 m 2 



a On^T(0,n) + a ln^T(2,n) + Ajn (^T(l,n) + ^T(l,n)) 



(A32) 



= +3M {£) 2 £ [(4o + #i* + 4s* 2 ) ST + (4o + d+t) (ST + SST)] ■ (A33) 

We have recently checked explicitly all the results of our 1997 paper and found out that they are 
equivalent with those of the present work if we make the approximation W = 2m. 

There are still two important sources of differences between these two sets of results. The first 
one is due to the fact that those of the earlier work were not truncated at a given order. The second 
one is that it did not include explicitly the two-loop diagrams, as we do now. In 1997 these effects 
were hidden within the irN sub-threshold coefficients and were therefore double counted. Even if 
these effects are numerically small, as we discussed in section IV of this work, this represents a 
rather important conceptual difference between both calculations. As two-loop contributions only 
arise at 0(q 4 ), the potential produced in 1997 would be numerically identical with the present one 
for distances larger than 1.5 fm if both of them were truncated at 0(q 3 ). 
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